Ultrafast demagnetization in the sp-d model: a theoretical study 
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We propose and analyze a theoretical model of ultrafast light-induced magnetization dynamics 
in systems of localized spins that are coupled to carriers' spins by sp-d exchange interaction. A 
prominent example of a class of materials falling into this category are ferromagnetic (III,Mn)V 
semiconductors, in which ultrafast demagnetization has been recently observed. In the proposed 
model light excitation heats up the population of carriers, taking it out of equilibrium with the 
localized spins. This triggers the process of energy and angular momentum exchange between the 
two spin systems, which lasts for the duration of the energy relaxation of the carriers. We derive 
the Master equation for the density matrix of a localized spin interacting with the hot carriers and 
couple it with a phenomenological treatment of the carrier dynamics. We develop a general theory 
within the sp-d model and we apply it to the ferromagnetic semiconductors, taking into account 
the valence band structure of these materials. We show that the fast spin relaxation of the carriers 
can sustain the flow of polarization between the localized and itinerant spins leading to significant 
demagnetization of the localized spin system, observed in (III,Mn)V materials. 

PACS numbers: 



I. INTRODUCTION 



During the last decade, ultrafast light-induced mag- 
netization dynamics has been a subject of many exper- 
imental studies, which have yielded results of interest 
from both the insight into the basic physical processes 
occurring in ferromagnets at short time-scales, and from 
the possible practical implications for ultrafast writing 
of magnetic memories. The seminal discovery^, of sub- 
picosecond demagnetization in Ni has challenged the pre- 
viously held belief that the fastest time-scale of mag- 
netization dynamics is the spin-lattice relaxation time, 
which is of the order of at least tens of picoseconds^. 
Since then, such an ultrafast magnetization quenching 
process in transition metals has been confirmed by var- 
ious time-resolved magneto-optical techinques 4 *^' 7 - as 
well al other experimental methods, such as magnetic sec- 
ond harmonic generation^ spin-resolved photoemission, 9 
and THz light emission 10 . Recently, an analogous phe- 
nomenon of light-induced demagnetization, including a 
complete destruction of ferromagnetic order in less than 
a picosecond, has been observed in (III,Mn)V ferromag- 
netic semiconductors i 1 1 1 12 

Apart from magnetization quenching, which is caused 
by absorption of light, there has recently been sig- 
nificant progress in nonthermal (coherent) light ma- 
nipulation of magnetic orderJ^ Coherent manipula- 
tion on the sub-picosecond time scale has been shown 
in ant iferromagneti o 14 ' 1 5 and ferrimagnetic materials^ 
Theories of coherent manipulation of magnetization 
by strong off-resonant light have been put forth 
for (III,Mn)V ferromagnetic semiconductors 17 and 
for closely related undoped paramagnetic (II,Mn)VI 
materials* 1 ^ In this paper we are interested in ultrafast 
demagnetization, which is an incoherent process involv- 
ing strong excitation. The laser pulse heats up the carrier 
population, and we analyze the magnetization dynamics 



induced by such a hot gas of electrons or holes. 

Ultrafast demagnetization experiments in transition 
metals were initially interpreted using a phenomenolog- 
ical three-temperature modeliiii 9 . In this approach the 
system is divided into three reservoirs: carriers, lattice, 
and spins. Excitation by light injects energy into the car- 
rier system. Each reservoir is described by its tempera- 
ture, and phenomenological equations are written down 
for heat flow between each pair of reservoirs. Physical 
underpinnings and time-scales of two of the couplings, 
carrier-lattice (electron-phonon) and spin-lattice interac- 
tions are known (for the latter see Ref. Q) . A new ingredi- 
ent, a direct carrier-spin coupling has to be postulated in 
order to explain the ultrafast demagnetization. In addi- 
tion to the lack of detailed microscopic understanding of 
such carrier-spin coupling, there are two major deficien- 
cies of the three-temperature model. The first is that 
treating carriers and spins as separate entities might not 
be a good starting point for the transition-metal itinerant 
ferromagnets. If such an approach is possible, the nature 
of an effective separation into subsystems should be elu- 
cidated. The second shortcoming is the fact that only the 
energy transfer between reservoirs is considered. As the 
key phenomenon to be explained is the change of mag- 
netization, a correct physical description should involve 
the mechanism by which the angular momentum (spin 
and orbital) is exchanged between the subsystems. 2 ^ 21 
Recently^ 2 Koopmans et. al. proposed a simplified model 
of a transition-metal ferromagnet in which spinless elec- 
trons induce spin-flips in a magnetic subsystem separate 
from the electronic system. Another proposed theory of 
demagnetization, 23 in which the cooperative effect of a 
coherent laser field and the spin-orbit coupling is calcu- 
lated, seems to be more suitable for the case of coherent 
manipulation of antifcrromagnetSi 13 i 24 

A model of ferromagnetism ideally suited for investi- 
gation of magnetization quenching induced by excitation 
of carriers is the sp-d model. In this approach, most 
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of the macroscopic magnetization comes from the local- 
ized d-shell spins (or f shells, in case of the rare earths), 
which are coupled by an exchange interaction to itin- 
erant s or p carriers. This model was introduced in- 
dependently by Zener— and by Vonsovskii 2 ^ in order 
to describe the transition metals, but was abandoned 
when the above-mentioned lack of sharp separation into 
s carriers and d spins was understood. Later, the s-d(f) 
model was succesfully applied to magnetic semiconduc- 
tors such as the chalcogenides of europium. 27 Recently, 
it has been again used in itinerant ferromagnets in order 
to analyze the situations in which an interplay between 
the transport and magnetic properties is critical. Ex- 
amples include spin-transfer torqu o 28 ' 29 and an enhance- 
ment of Gilbert damping in a magnet due to pumping of 
spin currents into adjacent non-magnetic material] 30 ! 31 
Finally, and most importantly for us, the p-d model is 
used to describe the basic physics of ferromagnetism in 
(III,Mn)V semiconductors ) 32 ! 33 such as Gai_ x Mn x As or 
Ini-^Mn^As . 

Ferromagnetic (III,Mn)V semiconductors are created 
by doping a small (~5%) molar fraction of Mn into a III- 
V semiconductor such as GaAs or InAs. The Mn ions, 
which substitute the cations, are acceptors and their d 
shells can be treated as well-localized 5/2 spinsi 3 ^ Hy- 
bridization between Mn d orbitals and anion p states 
leads to a kinetic p-d exchange interaction 3 ^., which cou- 
ples localized d spins to spins of p holes in the valence 
band. These holes are a source of indirect Mn-Mn spin 
coupling. At large hole densities (typically p«10 20 cm~ 3 ) 
a mean-field theory of Zene r 25 i 32 ' 36 has been successful 
in describing many features of ferromagnetism in these 
materials. It correctly predicts the critical temperature 
T c , increasing trend in T c with the density of carriers, 
and magnetic anisotropics. 33 

The strong correlation between the presence of the de- 
localized (or at least weakly localized) holes and ferro- 
magnetic ordering of Mn spins, has been firmly estab- 
lished experimentally. The optical induction of ferromag- 
netic transition (through cw photoinjection of holes) has 
been shown^ and the critical temperature and the co- 
ercive field have been altered by changing the density of 
holes in InMnAs by applying a gate voltage i 38 ! 39 Most 
relevant for this article is the observation of the sub- 
picosecond light induced demagnetization in InMnAs 11 
and GaMnAs^ A complete quenching of ferromagnetic 
order was achieved in InMnAs for pump fluences above 
10 mJ/cm 2 . This should be contrasted with behavior 
in Ni, where the demagnetization saturates at higher flu- 
ences while not reaching the complete demagnetization. 40 
Since the carrier concentration is much lower in ferromag- 
netic semiconductors as compared to metals, their mag- 
netization is more amenable to manipulation by external 
stimuli. 

These experiments led us to a preliminary investiga- 
tion of the spin dynamics induced by strong incoherent 
excitation within the framework of the sp-d model^ and 
to a proposal of the "inverse Overhauser effect" as the 



physical basis of the observed phenomenon. In this sce- 
nario, the rate of spin-flip scattering between the local- 
ized spins and the mobile carriers is enhanced by the non- 
equilibrium distribution of excited carriers. This triggers 
the transfer of angular momentum from the localized spin 
system to the carriers, leading to demagnetization of the 
localized spins. For this process to be sustained, the an- 
gular momentum transferred into the carriers' system has 
to be efficiently dissipated into the lattice by spin-orbit 
interaction assisted scattering, and we have stressed that 
the expected short spin relaxation time of holes in ferro- 
magnetic semiconductor is essential for the explanation 
of measured changes in magnetization. 

In this article, we present a detailed theory of ultrafast 
demagnetization within the sp-d model. We concentrate 
on ferromagnetic semiconductors, specifically, we discuss 
photoexcitation processes specific to these materials, and 
we perform calculations using an effective Hamiltonian 33 
model of the spin-split valence band. We assume that the 
indirect (carrier-mediated) exchange interaction between 
localized spins dominates over short-range antiferromag- 
netic d-d superexchange (it has been argued 3 ^ that in the 
presence of holes the superexchange is suppressed) . We 
also neglect the Mn interstitials 4 ^ which can form antifer- 
romagnetically coupled pairs with nearby substitutional 
Mn spins^ as we are only interested in Mn spins which 
participate in ferromagnetic order. Furthermore we work 
in the regime of strong carrier excitation, in which the 
carrier mediated Mn-Mn correlations beyond the mean- 
field level are averaged out. The localized spins flip inde- 
pendently in the common mean field due to the average 
carrier spin. We derive the rate equations for their dy- 
namics due to spin-flip scattering with the carriers spins. 
The transition rates depend on the instantaneous state 
of the carrier system at a given time. This derivation 
generalizes to the non-stationary case the theory from 
Ref. [43|, where the heating of Mn spins by electrons ex- 
cited by cw light was considered in a (II,Mn)VI based 
quantum well. Apart from the localized spin dynamics, 
we consider the energy and spin relaxation of carriers by 
means of phcnomcnological equations. The energy relax- 
ation time of the carriers sets an upper bound for the 
time-scale of the ultrafast demagnetization process. If 
the spin relaxation rate of carriers is smaller than the 
rate at which the angular momentum is injected into the 
carrier system by sp-d scattering, the demagnetization 
becomes suppressed. Below we show that this "spin bot- 
tleneck" effect does not affect qualitatively the results in 
p-type (III,Mn)V semiconductors. 

The general results of this article are directly ap- 
plicable to any system described by analogous s(p)- 
d(f) model. However, when the non-carrier-mediated 
d-d exchange coupling is strong (e.g. in the europium 
chalcogenides^), the starting point of the calculation 
should not be the interaction of a single spin with the car- 
riers (which is our approach here), but the enhancement 
of the carrier-magnon scattering. A similar treatment 
should be applied to a recent measurement^ 4 - of the ultra- 
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fast demagnetization in GaMnAs excited with the pump 
fluence 3 orders of magnitude smaller than in Ref. [ll|. 
In this case, our assumption of complete obliteration of 
the carrier-mediated exchange coupling by excitation of 
carriers might not hold, and spin-wave modes of coupled 
Mn spins should be the starting point of the calculation. 

The article is organized as follows: In section Hi Al we 
introduce the sp-d model. Section Hi Bl contains a discus- 
sion of the photoexcitation process, which is followed in 
Section Hl Cl bv the description of the model of the carrier 
bath heated up by absorption of light. Section |TTD] con- 
tains the formalism and the derivations of equations gov- 
erning the dynamics of localized spins interacting with 
hot itinerant spins. In Section HTT1 wc calculate the pro- 
cess of ultrafast demagnetization using a simplified band 
structure (a single spin-split band), introduce the equa- 
tions governing the carrier dynamics, and discuss how 
the energy and spin relaxation of carriers influences the 
demagnetization process. Complications introduced by 
valence band-structure of a ferromagnetic semiconduc- 
tor are discussed in Section HVi where we use a 6 band 
Luttinger Hamiltonian to calculate the hole-Mn spin-flip 
transition rate. Finally, we discuss the connection to the 
experiments in Section fVl 

II. GENERAL THEORY 
A. The sp-d model 

We use a Hamiltonian consisting of the single-particle 
carrier part He, localized spin part H$ and the sp-d 
coupling Hcs- 

H = Hc + Hs + H C s . (1) 

We do not explicitly consider the carrier-carrier interac- 
tion, carrier-phonon interaction, and the lattice dynam- 
ics. However, these interactions are included phenomeno- 
logically in our treatment of the carrier bath (Section 
III Cj) . The carrier part of the Hamiltonian is given by 

He = ^ erika^nkflnkcr , (2) 
nk 

where n is the band index, k is the wave- vector, and e nk is 
the energy. In the case of (III,Mn)V semiconductors, the 
valence band structure near the T point can be obtained 
from a k ■ p model, such as 6x6 Luttinger Hamiltonian 
commonly used in these systems.— Although the ferro- 
magnetic semiconductors are known to be heavily disor- 
dered, we have neglected the disorder potential in Eq. ([2]). 
If the disorder can be treated perturbatively, then its in- 
fluence on the treatment below is not expected to be 
strong. As we are working in the regime of very strong 
excitation of carriers, more subtle features of bands will 
be averaged out. On the other hand, if the effects of 
disorder are non-perturbative (e.g. the holes are in an 
impurity band rather than in the host valence band, as 



it was recently suggested in Ref. I45T ) , a new energy dis- 
persion (e.g. a single weakly dispersive band) should be 
included introduced in Eq. J2]). 

The localized spin Hamiltonian Hg includes Zeeman 
coupling to the external magnetic field. In the ferro- 
magnetic phase, the splitting of the localized spin S due 
to exchange interaction with carriers exceeds the typical 
Zeeman splitting, and the only role of Zeeman term is to 
choose the direction of the magnetization. We disregard 
possible S-S exchange coupling by mechanisms other 35 
than indirect carrier-mediated exchange. 

Finally, the exchange coupling is given by 

Hcs = ~Yl E S ; >k|s>'kV( k '- k ) R < at„ k a n , k , , 

(3) 

where S/ is the spin operator of localized spin at R; , § is 
the carrier spin operator and 7 is the exchange constant. 
In the literature on diluted magnetic semiconductors, the 
7 parameter is called a and (3 for the conduction and 
valence band electrons, respectively^ The typical values 
in (III,Mn)V materials are awlO and /3~— 50 meV-nm 3 . 
The exchange energy J per unit cell is J=N j, where N 
is the density of the cations. Accordingly, Aoa«0.2 eV 
and N f3&-1 eV. 

For a single band, the sum over n states is simplified 

to 

^cs = -^EE^ r §kk ' e 4(k '~ k)Ri , (4) 

l k,k' 

where the scalar product of spin operators can be ex- 
pressed using spin ladder operators as: 

S, ■ s kk , = i(£+s kk , + Sfs+v) + SfS kk , , (5) 
and the explicit form of carrier spin operators is 

*kk' = o( a W«k'+ - e^k-ak'-) (6) 



«kk' = ftt k+ a k'- , s kk , = a t k _a k / + . (7) 

The magnetization (average localized spin) is in the z 
direction, and we treat the sp-d Hamiltonian in a mean- 
field approximation. The mean-field felt by the carriers is 
obtained by performing the simplest virtual crystal dis- 
order averaging of localized spins positions: 

H™ f = -nij(S z ) 2J (?ik|s z |n'k)a t nk a„/ k , (8) 

n,n',k 

where Ui is the density of localized spins and (S z ) is the 
average localized spin. The energy of the carrier spin 
splitting is defined as 

A = -nMS z ) . (9) 
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The typical value of A in GaMnAs and InMnAs with the 
highest critical temperatures is about 0.1 eV. By averag- 
ing Eq. ([3]) with respect to the carriers we arrive at 

Hf=Sj2Sf, (10) 

i 

where the energy splitting of different m levels of a local- 
ized spin is 

5 = -n cl (s z ), (11) 

with n c the carrier density and (s z ) the average carrier 
spin. The typical magnitude of S in (III,Mn)V materials 
is of the order of 1 meV. In the mean-field theory of 
ferromagnetis m 32 ' 36 ' 47 only the two averaged fields (S z ) 
and (s z ) are taken into account, and the free energy of 
their interaction is minimized. 



B. Light excitation of the carrier system in 
(III,Mn)V 

To construct a model of the excited carrier bath, it is 
helpful to first analyze qualitatively the process of carrier 
photoexcitation in ferromagnetic semiconductors. The 
cw (magneto)optical spectra of GaMnAs are qualitatively 
different from the ones in pure GaAs, as they show very 
strong effects of disorderi 48 i 49 i 50 i 51 In GaMnAs, there is 
no gap in the absorption, which is at least of the order of 
10 4 cm -1 for all energies^ The origin of strong absorp- 
tion inside the host gap is a matter of controversy (see 
Figure [1]). It has been proposed that the additional opti- 
cal transitions are between the valence band and the dis- 
persionless levels (located about 0.7 E g above the valence 
band) originating from As antisites or Mn interstitials. 52 
Transitions terminating inside the Mn derived impurity 
band ~0.1 eV above the valence band edge have also been 
suggestedi^^^ 3 - 

In addition, the role of the inter-valence band tran- 
sitions is increased in the disordered material, in which 
the k-selection rule is relaxed. As a result, the absorp- 
tion of light with energy smaller than the bandgap of the 
host material leads to a strong excitation in the valence 
band, which occurs due to inter-valence transitions and 
possible transitions from below the Fermi energy into the 
localized states within the gap. The initial distribution 
of photoholes is expected to be very broad, determined 
only by energy conservation (not by k-selection). We use 
the broadness of the distribution of the carriers after the 
excitation and their total number as free parameters of 
the theory. Some flexibility in the used values of total 
hole density immediately after the excitation p is also 
justified by the fact that the initial densities po are usu- 
ally known with only an order of magnitude accuracy. 
We use an energy width of hole distribution after exci- 
tation of the order of 0.1 eV, a value comparable to the 
exchange splitting of the valence band. In Section [V] we 
use the experimental data to argue that the above energy 
scale is sensible. 
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FIG. 1: (Color online) The band structure of GaMnAs (solid 
lines) calculated using 8 band k • p model with mean-field 
sp-d exchange interaction (|A|=0.15 eV). The dotted lines de- 
note the possible positions of midgap energy levels of different 
origins. In the model without an impurity band, the Fermi 
energy is typically around —0.1 eV for hole density p~10 20 
cm -3 . 

The optical experiments in InMnAs are limited to sam- 
ples with very small Mn concentrations— or energies 
above the fundamental bandgap^ The role of various de- 
fects has not been investigated in InMnAs as closely as in 
GaMnAs, but the observed^ ultrafast electron trapping 
time suggests that InMnAs does have a large concentra- 
tion of midgap defects characteristic of low-temperature 
grown III-V semiconductors. We can also use the results 
of this article to argue for the presence of strong transi- 
tions involving the valence band states and states not in 
the conduction band, on the ground that the number of 
states available in the conduction band is far too small 
to explain the demagnetization results. 

This picture of photoexcitation is confirmed by the ul- 
trafast demagnetization measurements^ in GaMnAs ex- 
cited by 0.6 eV pump, far below 1.5 eV bandgap of GaAs . 
The results are similar to the ones observed^ in InMnAs 
(bandgap of 0.4 eV) excited by the same pump beam, 
showing that excitation into the conduction band plays 
a minor role in the demagnetization process. 



C. Model of the carrier bath 

According to the discussion in the previous section, 
immediately after the photoexcitation the distribution 
of carriers is very broad. We will approximate it by 
a thermal distribution described by the electronic (car- 
rier) temperature T c . Using the experimental data from 
Ref. [TT| we estimate the initial T c to be of the order of 
1000 K (for details see Section EJ. 

The interaction between the carriers and the localized 
spins produces a mean-field term (Eq. ((5J)) and a sec- 
ondary term corresponding to simultaneous flips of the 
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itinerant and localized spins. The latter causes exchange 
of angular momentum between the excited carriers and 
the localized spin system. Due to the spherical symme- 
try of the sp-d Hamiltonian, the sp-d interaction alone 
conserves the total spin, and it can only move the spin 
polarization from one system to another. The mecha- 
nism of this transfer of spin is described in detail in the 
next Section. Now we concentrate on the features of the 
carrier bath which are specific to the case at hand: the 
possibility of dynamic spin polarization and the energy 
relaxation of carriers. 

The spin transferred into the carrier system is not con- 
served due to the spin-orbit interaction. In its presence 
the scattering within the electronic system is accompa- 
nied by spin relaxation 56 (ultimately into the lattice). 
However, the spin relaxation occurs on a finite time-scale 
r sr . If the rate at which the carriers relax their spins is 
smaller than the rate at which angular momentum is in- 
jected into the carrier population, there is a dynamic spin 
polarization of the carriers. In this case the average car- 
rier spin deviates from the mean-field value determined 
by A splitting of the bands. In the case when the car- 
riers occupy a single spin-split band, there is a simple 
way of introducing the dynamic polarization into our for- 
malism (the case of multiple bands strongly coupled by 
spin-orbit interaction is considered in Section IIV[) . We 
assume that populations for both spin directions are de- 
scribed by Fermi-Dirac distributions with the common 
temperature T e but different quasi-chemical potentials 
/j, s (spin s=|,|). This is equivalent to saying that we 
coarse-grain the dynamics on a time-scale larger than the 
time in which the energy is redistributed within the car- 
rier system (by carrier-carrier scattering) , and we explic- 
itly consider only the slower processes: spin relaxation of 
carriers and energy transfer into the lattice. It is simi- 
lar to the situation encountered in semiconductor lasers, 
where the processes of thermalization of electrons and 
holes separately occur faster than the recombination, and 
the resulting globally out-of-equilibrium situation can be 
described by using the quasi-equilibrium form of the car- 
rier distribution, with different quasi-chemical potentials 
for electrons and holes. 

As we discuss in detail in the following sections, dy- 
namic polarization of carriers affects the spin-flip scat- 
tering rate. This is a back-action effect of the S system 
on the carrier system: demagnetization of localized spins 
causes the dynamic polarization of carriers, which in turn 
influences the rate of transfer of angular momentum be- 
tween the systems. Additionally, as the average S spin 
changes in time, so does the A splitting of the carrier 
bands. 

Apart from the dynamic polarization effects, we also 
have to take into account that the carriers are not in 
thermal cqulibrium with the lattice, and in the first pi- 
coseconds after the photoexcitation they are described 
by a different temperature than the lattice. We model 
the carrier-phonon interaction (leading to the cooling of 
carriers) and carrier-carrier scattering (maintaining the 



thermal distribution) phenomenologically. In metals, the 
electron-phonon energy relaxation time is a couple hun- 
dreds of femtoseconds , 57 i 58 In semiconductors, the regime 
of such strong excitation as used in Ref. |Tl| to demagne- 
tize InMnAs has not been investigated in detail. The- 
oretical calculations of energy relaxation of holes deep 
in the valence band indicate that emission of optical 
phonons is very efficient for these states^ Calculated 
emission of 7 phonons during a picosecond should cor- 
respond to significant cooling of the carriers. Naive ex- 
trapolation of energy-loss curves^ to hole temperatures 
of ~1 000K also gives a sub-picosecond energy relaxation 
time. In the following, we will use an energy relax- 
ation time te of the order of picosecond, and assume 
T c (t) — T c (0) exp(— ^/te). More generally, the carrier 
and lattice temperatures could be solved for using a two- 
temperature model.— However, in the case of (III,Mn)V 
semiconductors the rise in lattice temperature due to 
heat transfer from the carriers described by such T c is 
quite small, so that the final common temperature of 
carriers and the lattice is less than 100 K. As we discuss 
in detail further on, such a temperature corresponds to 
very slow demagnetization dynamics, and for the purpose 
of calculating the ultrafast process, we can then assume 
that the carrier temperature simply decays towards zero. 

In the following, all the averages with respect to carrier 
degrees of freedom will be taken using a density matrix 
of non- interacting electrons pc , which is not necessarily 
of the equilibrium form, but remains diagonal in the ba- 
sis of carrier Hamiltonian's eigenstates |a)=|nk). These 
averages are denoted as (...)c=Ttlc{pc---}i and they are 
given by 

(a t Q a /3 )c = S a p f a (12) 
(cr ' aapa) \as)c = 5 a p8 lS / Q / 7 (13) 

+S a sS/3 7 f a (l - ffj) 

where f a is the average occupation of a state. In a single- 
band model, the a.. .6 indices refer to |sk) states with spin 
,s=|,j. /sk is a Fermi-Dirac function at temperature T e 
with spin-dependent chemical potential /i s . 

D. Rate equations for the localized spin 

The framework of the sp-d model allows for clear sepa- 
ration of carrier (C) and localized spin (S) systems. The 
mean-field parts of their mutual interaction are given by 
Eq. {SJ and (jTUJ) . The spin-splitting A of the carri- 
ers' band is proportional to the instantaneous (S z ), and 
the splitting of localized spins, 5, changes with the av- 
erage carrier spin (s z ). We assume that any correlation 
between localized spins beyond the mean-field Zener ap- 
proach is obliterated by the strong scattering of excited 
carriers. Each localized spin feels the dynamics of the 
other spins only through their average value, which in- 
fluences the state of the carrier system (its spin splitting 
A). Essentially, we consider an ensemble of paramag- 
netic S spins interacting with a bath, the properties of 
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which depend on the average S. The ferromagnetism en- 
ters only as an initial condition: the S system is polarized 
at t=0. Below we derive the equations for the dynam- 
ics of the average S spin due to the interaction with the 
carrier bath, the state of which depends on the average 
S. 

The Hamiltonian of the localized S spins and the car- 
riers system is written as 

H s -c = H + V = H C + H™ f + ]T(<5Sf + V t ) , (14) 

i 

where He is the carrier band Hamiltonian (Eq. ©), 
H™f is the mean-field spin-splitting from Eq. {8]), 6 is 
the mean- field localized spin splitting defined in Eq. (fTTjl , 
and the spin-flip term of l-th localized spin V\ comes from 
part of sp-d interaction which is off-diagonal in S z basis. 
We write it in the following way: 

Vi = S+F- + S^F+ (15) 

where F ± are proportional to the ladder operators of 
carrier spin. In the general case of multiple bands (as in 
Eq. ©) we have 

P± = ^ (nklsVkVnkan'k' , (16) 

nk,n'k' 

whereas for a single band we have 

F ± = -^T, a ^'T ■ ( 17 ) 

k,k' 

Now we follow a standard way of deriving the Master 
equation for the density matrix of the localized spin sys- 
tem interacting with a carrier bath (see e.g. Ref.l62l). The 
total density matrix of the system is assumed to factorize 
into the carrier and the localized spin density operators: 

p(t) « pc{t)p 8 {t) , (18) 

and the Liouville equation for time-dependence of p(t) 
is turned into an equation for ps(t) by tracing out the 
carrier degrees of freedom. However, unlike in standard 
treatment^ the state of the carrier bath changes in time, 
as discussed in the previous section. 

The usual derivation of the Master equation implies 
coarse-graining of the system dynamics on time-scale At 
longer than the correlation time of the bath, r c . This is 
the condition on which the Markov approximation rests. 
For the gas of carriers described by the effective tem- 
perature fcsT e sa0.1 eV, we can expect this time to be of 
the order of a few femtoseconds. Let us define a time- 
scale r p in which changes in the carriers' density matrix 
pc(t) occur. The contributions to the evolution of pc(t) 
are as follows. The carrier temperature changes appre- 
ciably during energy relaxation time te~1 ps. The spin 
splitting of the band A(t) is proportional to the average 
localized spin (S(t)), which we expect to decrease during 



a characteristic demagnetization time tm- The build-up 
of the dynamic spin polarization of the carriers is de- 
termined by two processes: transfer of spin from the S 
system occurring during the aforementioned time tm and 
the spin relaxation of carriers characterized by time r sr 
(for very short T sr , there is no dynamical polarization). 
Now, we will assume that all these result in the time r p 
much larger than the bath correlation time r c , so that 
we can choose our coarse-graining step A fulfilling the 
condition: 

t c < At < t p . (19) 

In such a case, at each coarse-grained time-step t n we 
can derive a Master equation with carriers described by 
Pc{tn) treated as approximately constant during At. In 
this way we can use the Markov approximation locally 
in time, having separated the "macroscopic" back-action 
of the localized spin system on the bath, which occurs 
on a longer time-scale. The Master equation for the lo- 
calized spin density matrix 5q is then derived exactly as 
in the usual case (see Ref. I62I ). only with the transition 
rates depending on time t n . We write the equations in 
the continuum limit, keeping in mind that they cannot 
be used for times shorter than At, thus obtaining the 
following rate equations for the diagonal elements of the 
localized spin density matrix p^=p^ m : 

J t Pm = -(W m -l,m + W m +l,m)Pm 

+w m , m +ipi+i + w m , m -ipi-i , (20) 

where W n<m is the transition rate from m to n energy 
level of the localized spin induced by sp-d interaction 
with the carriers. The time dependence of the transition 
rates is suppressed for clarity. Then the average localized 
spin evolves according to 

!«.(*» -5>i<4. < 21 > 

m 

The general formula for W m ^ m ±i is 

1 f°° 
W m , m±1 {t) = -S^ m±1 J^ dt'e ±M '^C T± (t;t') , 

in which the matrix element (squared) for the flip of the 
localized spin is given by 

S^ m±1 = S(S+l)-m(m±l) , (23) 

where S is the magnitude of the localized spin. The cor- 
relation function of the carriers Cy(t;t'), with i,j=±, is 
given by (reverting to the coarse-grained notation with 
t n replacing t) 

Cy(i n ;t') = Tr c {pc{tn)K{t n + t')Fi{t n )} , (24) 

where Trc{...} is the trace with respect to the carrier 
states and F^(t) are operators defined in Eqs. (TTS|) and 
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(fT7| written in the "local" interaction picture at coarse- 
grained time t n : 

Fi(t") = exp{~#o(in)*"} F l (t")exp j-^o^)*" 

(25) 

where F l {t) is the operator in the Schrodringer picture, 
and Ho(t n ), defined in Eq. (fH|) , depends on time through 
the mean- field spin splittings S and A. These correlation 



functions Cij(t n ;t') decay for t' larger than the correla- 
tion time t c , which has to be much shorter than the time 
on which pcif) changes. For this reason the integration 
domain in Eq. (f2"2")l is effectively t' € (— t c ,t c ). In this 
range of t', the above definitions of Cy and F^(t) make 
sense. 

In the general case of the multiple bands the transition 
rates are given by: 



W„ 



±1=JY S ™MiE| (^3 / 7^3l( n ' k 'l J± l nk )| 2 /nk(l - fn'k')S(ink ~ £„'k' ±5) , 



(26) 



r 



where n and n' are labeling the subbands, /„k is the oc- 
cupation of nk) state, and e n k are the band energies with 
the mean-field exchange interaction with localized spins 
taken into account. The distribution functions, energies 
e„k, and S depend implicitly on time, as we discussed 
before. 

For a single spin-split band we replace n and n' by two 
spin indices and recover the formula given in Ref. [ill. 
An analogous expression has been used in Ref. l43l where 
heating of the Mn spins by photoelectrons was consid- 
ered in a paramagnetic (II,Mn)VI quantum well. Actu- 
ally, in the case of carriers being a true reservoir of energy 
and polarization, Eqs. (|2"0|) and (j2"B"|) were derived orig- 
inally by Korringa 5 - in order to describe the relaxation 
of nuclear spins coupled to carriers' spins by hypcrfinc 
interaction. 

The rate equation (|2H)) has to be augmented by equa- 
tions governing the dynamics of the carrier distribution 
function /„k- A discussion of these additional equations 
follows in the next section. 



III. DEMAGNETIZATION DUE TO CARRIERS 
IN A SINGLE SPIN-SPLIT BAND 



Let us concentrate now on a model of a single spin- 
split band of s symmetry. From it, we deduce the general 
features of the behavior of the system in a simple way. 
The treatment of the dynamical polarization of carriers 
is especially transparent in this case. We make use of 
the distribution functions f s for spin s=|,J. characterized 
by two different chemical potentials and a common 
temperature T e (see Sec. Ill C|) . 

We define reduced transition rates W+- and W- + 
given by W m ^ m -\ and W m>m +x, respectively, with the 
localized spin matrix elements 5„ m=F i removed (see 
Eq. (|2"6"|) ). In a single band with spin-splitting A, the 
transition rate W-\ can be rewritten using the spin- 



resolved densities of states D S (E) 
7 2 2tt 



W+- 



4 h 



xD 1 (E)D l (E-S) 



(27) 



and W—+ is obtained by exchanging the spins and chang- 
ing the sign of S. Manipulating the explicit forms of oc- 
cupation functions we obtain a generalization of the de- 
tailed balance condition to the case of different chemical 
potentials for two spin directions: 



W. 



W+- 



(28) 



where /3 =l/fcsT e and A/z=//f — (Xi is the spin splitting 
of the carriers' chemical potential. For /i|=/i| we recover 
the usual detailed balance condition. 

When 5=1/2, Eq. (f2"Tj) can be transformed into the 
Bloch-like equation for the dynamics of average localized 
spin (S z (t)): 



d , Q( ^_ (S z (t)) - S (t) 



(29) 



where 5o(t) is the instantaneous equilibrium value of the 
spin, given by the transitions rates at time t: 



^0 = 77 i 



21^4 



-^± = -itanh(/3 e (5-A/x)) , (30) 



and the relaxation time is given by 

T 1 (t) = (W+_+W_ + )- 1 



(31) 



Note that in applications to (III,Mn)V semiconductors 
we are going to be interested in the regime of /3 <5<§;1 and 
in the localized spin 5=5/2. If the dynamic spin splitting 
also fulfills /3 c A/i<g;l, we can approximate Eq. ((21)) for 
any magnitude of spin S by the following expression: 



dt 



(S z (t)}~-2W+^(t)(S z (t)) 



(32) 
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For the conditions considered below, this equation gives 
a very good description of the initial stage of the lo- 
calized spin dynamics, in which the carrier temperature 
is very high. When the temperature drops so that the 
above inequalities are violated, one has to solve the full 
eqs. ([20]) and (J2TJ) for S>l/2, and the Bloch equation 
for S=l/2. 

We now introduce the phenomenological equations de- 
scribing the dynamics of the carrier bath. The time de- 
pendence of T e comes from the cooling of carriers by 
phonon emission, and we model it by a simple decay (see 
Section [EC]) 

T c (t) =T c (0)e-*/ rE , (33) 

where te is the energy relaxation time of highly excited 
carriers. The changes of chemical potentials are governed 
by a second phenomenological equation for the dynamics 
of the average carrier spin (s(t)}: 

d {s{t)) = _^d _(s(t))- So (A,T c ) 
dt n c dt T sr 

The first term on the right describes the transfer of an- 
gular momentum by spin-flips (with n c /rii being the ra- 
tio of the carrier density to the localized spin density), 
and the second term describes the relaxation (on time 
scale r sr ) of the average carrier spin towards the instanta- 
neous equilibrium value sq determined by the spin split- 
ting A and the carrier temperature T c . The above set 
of single-band equations was used in Refs. [ll| and [H 
to qualitatively model the demagnetization in (III,Mn)V 
semiconductors j 11 ' 12 



A. The inverse Overhauser effect 

The demagnetization process described by the above 
equations occurs in the following way. We model the ab- 
sorption of a ^100 fs light pulse as an instantaneous in- 
crease of concentration and temperature of carriers. The 
heating of the carriers by a pulse of light modifies the 

spin- flip transition rates W-\ and W (.. The broader 

the carrier distributions are (the higher the T e is), the 
larger these rates are. When the dynamic spin splitting 
A/i is zero, the detailed balance of transition rates tells 
us that the localized spin is going to evolve towards a 
new value corresponding to a high temperature T e . This 
final value of (S z ), and the rate at which it is approached, 
change with the decrease of the carrier temperature. In 
addition, the possible build-up of polarization of the car- 
riers changes the spin-flip transition rates in such a way 
that the spin transfer is blocked, and without spin relax- 
ation in the carrier system a "polarization bottleneck" 
occurs. This is analogous to the "magnetic resonance 
bottleneck" known in electron spin resonance of local- 
ized moments in metals.— In the latter, the resonance 
spectrum of the localized spins is changed when the spin 
relaxation of carriers is long enough for the two spin sys- 
tems to become "locked" together in precession. Here, in 



the extreme case of very slow spin relaxation, the initial 
demagnetization can result in a flip of all the carriers of 
one spin direction, leaving the carrier system in a "dy- 
namic half- metallic state" . After such a saturation of de- 
magnetization process, the rate of spin flip is determined 
by the (low) rate of the carrier spin relaxation. 

The basic principle of demagnetization is analogous to 
the well-known Overhauser effect^ in which the itin- 
erant spins are optically pumped, and this injected po- 
larization is transferred by s-d type interaction to the 
localized spins. Although the original Overhauser ef- 
fect involves pumping of angular momentum into one of 
the spin populations, the essence of the effect is taking 
one spin population out of equilibrium with another and 
thus inducing the transfer of angular momentum between 
them. One generalization of the Overhauser effect was 
proposed and realized experimentally by Fcher. 66 i 67 The 
idea was to heat up the electrons by passing a current 
through the sample in order to induce spin-flips between 
the electrons spins and nuclear spins. Depending on the 
parameters of both spin systems, increase of either nu- 
clear spin polarization 6 ^ or electron polarization 68 is pre- 
dicted. The "inverse Overhauser effect" presented here 
is related to the latter: by heating up the carriers, we 
induce the transfer of angular momentum from localized 
spins to electron (hole) spins. 



B. Carriers as a reservoir of angular momentum: 

We first analyze the case in which the carrier spin re- 
laxation is so fast that the carriers are a good sink of po- 
larization, so that A[i=0 at each moment of time. The 
occupation factors f s are the same for both spins, and 
are characterized by the time varying temperature T c and 
chemical potential fi. 

When the temperature is not too high, i.e., fcgTe is 
smaller than the energy scale on which the densities of 
states change appreciably (but still much larger than lo- 
calized spin splitting 5), the transition rate (|2"T|) can be 
approximated by 

W+-(t) « 2_^fe B T e (i)£> + (jm) . (35) 

It shows that the rate of demagnetization scales with T e 
and that large densities of states for both spins around 
the Fermi level are needed. Together with the 7 2 scaling, 
this shows that in (III,Mn)V ferromagnetic semiconduc- 
tors the holes (having larger mass and exchange constant) 
are much more effective in the demagnetization process 
than the electrons. 

Fig. [2] and [3] illustrate the dependence of W-\ from 

Eq. (|27p on temperature T c , spin splitting A, and carrier 
concentration. We use the density of states of a parabolic 
band, with effective mass m ff=l, which roughly corre- 
sponds to the density-of-states mass in the valence band 
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FIG. 2: (Color online) Transition rate in a sin- 
gle parabolic band with effective mass m c g—l and 7=50 
nieV-nm 3 as a function of carrier temperature and band split- 
ting A. The concentration of carriers is 10 20 cm -3 . 
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FIG. 3: (Color online) Transition rate W+- for A=0.15 eV 
and different carrier temperatures, as a function of carrier 
density. Other parameters are the same as in Figure [5] 



of GaMnAs. W_| goes to zero for low carrier temper- 
ature, when phase-space blocking limits the number of 
states which can scatter. It also decreases for increas- 
ing spin-splitting A, when the number of minority spins 
available for spin-flip goes down. This effect is stronger 
for smaller carrier concentrations and lower tempera- 
tures. For T e close to 1000 K, the corresponding charac- 
teristic time Xi is of the order of a picosecond. Whether 
a full demagnetization occurs depends on the rate of the 
energy relaxation of carriers. If the carrier temperature 
does not drop significantly within the time Ti(0) (calcu- 
lated at the initial carrier temperature T c ), then a sig- 
nificant demagnetization occurs on this time-scale. On 
the other hand, if the T e changes strongly on the scale of 

Ti(0), then we have to solve our equations with (t) 

updated according to carrier temperature changes from 
Eq. (|33|) . We are interested in the dynamics occurring 
during the first picosecond. When T\(t) becomes much 
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FIG. 4: (Color online) Demagnetization of the localized spins 
S=5/2 (normalized to 1) in a single band model with m e ff=l 
and 7=50 meV-nm 3 for different energy relaxation times of 
carriers. The concentration of carriers n c =10 20 cm -3 , and 
the density of localized spins is such that the band split- 
ting A=150 meV (corresponding to i«0.055 in ferromagnetic 
semiconductors) . 



larger than 1 ps, then for our purposes the demagnetiza- 
tion process is effectively stopped. Such an effect of cool- 
ing of carriers on demagnetization is illustrated in Figure 
31 For all energy relaxation times te the initial slope 
of (S z (t)) is the same, given approximately by Eq. ([55]) 
evaluated for T c (0), but the time at which the demagne- 
tization ceases and the saturation value depend on te- 
This shows that the time-scale on which the ultrafast de- 
magnetization occurs can be given by te, which is not 
related to magnetic properties of the material. From this 
point of view, it is not the fact that the magnetization 
drop occurs in less than a picosecond which is interest- 
ing. Instead it is the magnitude of the demagnetization 
which demands explanation. 



C. The effects of the dynamic spin polarization of 
carriers: A/i^O 

The efficiency of the "inverse Overhauser effect" is lim- 
ited by the finite spin relaxation time of the carriers. 
Each spin flip which leads to the demagnetization, if not 
followed by carrier spin relaxations, diminishes the phase 
space available for next transition of this kind. The re- 
sult is the decrease of the net number of spin-flips during 
the time te, which translates into smaller total demag- 
netization. 

Let us concentrate on the case in which W-\ is the 

transition which leads to the demagnetization of the lo- 
calized spins. The corresponding electron spin flip is from 
spin | to |, so that the demagnetization of S spins leads 
to A/i<0. In the limit of negligible S and A/j, smaller 
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FIG. 5: (Color online) Demagnetization of the localized spin 
S for different carrier spin relaxation times, with energy relax- 
ation time tb=0.5 ps. The effective mass of a single spin-split 
band is a) m e ff= 1 and b) m c ff=0.5. All the other parame- 
ters are same as in Fig. [4] The total drop in magnetization 
becomes smaller for longer carrier spin relaxation times. 



than ksT c , we can approximate Eq. (|35|) by 



W+-(t) 



~/ 2 2ir 

~ fc S T e (t)£> + (Mi)^-(M|) 



X 1 + 



A^ 

2k B T c 



(36) 



From this we see that W-\ decreases in the presence of 

dynamic spin polarization. When |A/z| becomes compa- 
rable to ksTc, W.\ goes to zero. 

A rough estimate of what r sr is short enough to be con- 
sidered instantaneous can be given assuming a constant 
density of states D for both spins. The localized spin 
splitting is then |<5| = |7A|_D/2, and the second term in 
Eq. (|34[) . corresponding to spin relaxation can be writ- 
ten as DA/i/2n c T sr . We want it to dominate over the 
first term (spin transfer from localized S) for A/z small 
enough, that the transition rate W+_ is still unaffected 
by such dynamic spin splitting, i.e. for Afj,^ksT c . The 
resulting inequality is 



Dk B T c 

iSm 



>W+-(Afx^0,T e )r sl 



which means that the ratio of the density of carriers avail- 
able for the spin flip to the localized spin density is larger 
than the product of spin-flip rate and carrier spin relax- 
ation time. Using Eq. (|35[) we can transform this in- 
equality into a very simple, but physically less intuitive 
form: 



T sr < 



h 

47T<5 



(38) 



In (ni,Mn)V semiconductors for typical value of local- 
ized spin splitting 5 is of the order of 1 meV, the spin 



relaxation time has then to be smaller than 100 fs for 
spin bottleneck to become unimportant. 

The "spin bottleneck" effect is illustrated in Figure 03 
where calculated (S z (t)) are plotted for different values 
of carrier spin relaxation time r sr . The bottleneck effect 
is stronger in Figure [5^, where the effective mass m e s=l 
and carrier concentration n c =10 20 cm~ 3 . In Figure [SJa, 
where m e ff=0.5, the difference between results for r sr =10 
fs and r sr =0 is smaller. When the effective mass m e g is 
smaller, while keeping n c and A the same, the average 
carrier spin is decreased, as the Fermi energy becomes 
larger compared with A. Then 5 is smaller, making the 
inequality (|38|) easier to fulfill. 



IV. 



DEMAGNETIZATION DUE TO HOLES IN 
THE VALENCE BAND OF (III,MN)V 
SEMICONDUCTOR 



For the case of (III,Mn)V semiconductors the carriers 
relevant for ultrafast demagnetization are the holes. If 
they reside in an impurity band, the single band theory 
described above could be applied, but currently there is 
no simple quantitative model for this case. Below we 
perform calculations for the case of holes residing inside 
the valence band of a bulk semiconductor. We use an 
"effective Hamiltonian" mode l 33 ' 36 in which a mean-field 
p-d term from Eq. © is added to the 6 band Luttinger 
Hamiltonian. In this way we can analyze quantitatively 
the influence of the strong spin-orbit interaction on the 
spin-flip transition rate. 

A typical plot of energy dispersions in (III,Mn)V cal- 
culated by k • p method with exchange splitting A=0.15 
eV is shown in Figure[T] The corresponding spin-resolved 
densities of states are shown in Figure [6l In order to cal- 
culate the W + _ transition rate we cannot use Eq. (|27[) . 
which requires only the densities of states. In the p sym- 
metry band the matrix elements of carrier ladder oper- 
ators are non-trivial due to spin-orbit coupling, and 
cannot be neglected. The full Eq. (|2"6"|) has to be used. 

Strong spin-orbit interaction suppresses the spin-flip 
transition rate, compared to the value predicted by 
Eq. (|27|) . If we disregard for a moment the orbital 
parts of the |nk) eigenstates, the squared matrix element 
|(n'k'|s~|?ik)| 2 is maximal when the spin part of |nk) 
(37) state is If), and the spin part of |n'k') is | |). If both 
states contained equal mixtures of spin up and down, the 
squared matrix element would be 1/4. When the orbital 
parts of the states are taken into account, the spin-orbit 
interaction diminishes the matrix elements between the 
states with different spin character, because it aligns or- 
thogonal orbital wave functions (of different orbital an- 
gular momentum I) with opposite spins s. In the limit 
of infinite spin-orbit interaction, a spin-flip between pure 
spin-up and spin-down states is impossible, because the 
orbital parts of |nk) states are exactly orthogonal. This 
limit is realized in the case of the 4 band Luttinger Hamil- 
tonian for the holes. We have evaluated analytically the 
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FIG. 6: (Color online) Valence band spin-resolved densities 
of states in GaMnAs. The solid lines are for single spin with 
A=0 in GaMnAs and InMnAs. Dashed and dotted lines are 
majority and minority spin for splitting A=0.15 eV in GaM- 
nAs. 



(a) 




1 1.5 2 2.5 3 



hole density (cm 3 ) xio 20 

FIG. 7: (Color online) (a) Reduced transition rate W+- in 
GaMnAs for hole density p=10 20 and no dynamic polariza- 
tion of holes (Ap=0). Analogous plots for other values of p 
and for InMnAs look qualitatively similar, only the overall 
magnitude of W+- changes, (b) Hole density dependence of 
W+- at T e =1000 K and spin splitting A=-0.15 eV (the lines 
are guides for the eye). 



spin-flip transition rate in a spherical 4x4 Luttinger 
model at zero spin splitting, and obtained that the result 
of the exact Eq. (|26|) is smaller by a factor of 5/18 than 
the value obtained from densities of states in Eq. (|27[) . 
This sets a limit on how much the spin-orbit interaction 
can suppress the transition rates for small A in a full 6 
band model. 

For the actual calculation of the time-dependence of 
magnetization, we have evaluated the transition rates 
from Eq. (j2"6|) using the band-structure obtained from the 
6x6 Luttinger Hamiltonian. In all the following calcula- 
tions, we have put 5=0, which is justified by the small- 
ness of (5 in ferromagnetic semiconductors and the fact 
that we are primarily interested in the regime of high 
T c , where 6 can be completely disregarded. In Figure [7| 
we plot W-| for GaMnAs and InMnAs for various car- 
rier temperatures, exchange splittings and hole densities. 
The transition rates are larger for GaMnAs, which can 
be traced mostly to larger density of states in the valence 
band compared to InMnAs. 

The presence of strong spin-orbit interaction makes in- 
cluding the dynamic spin polarization of carriers much 
harder than in a single-band case. Due to a large spin- 
orbit splitting Aso— 0-3 — 0.4 eV the |nk) states have 
mixed spin character even in presence of exchange split- 
ting A«0.1 eV, and a simple introduction of different 
distribution functions f s for two spin directions is not 
possible. Performing an ensemble Monte Carlo calcula- 
tion of dynamics of strongly excited holes^ with spin de- 
gree of freedom taken into account is a vast undertaking 
beyond the purpose of this investigation. We resort to an 
approximate method of including the effects of non-zero 
Afi. 

Above we have calculated the exact transition rates 
for A/i=0. We have also calculated the corresponding 
W-f assuming constant matrix elements and only tak- 
ing spin-resolved densities of states calculated from Lut- 
tinger Hamiltonian (Eq. ([27]) ). The inclusion of matrix 

elements leads to decrease of by a factor of 0.25-0.3 

in the most relevant range of high T c and A~0.1 eV (not 
shown). We assume that this ratio of exact and DOS- 
derived W-| also holds in the case of a finite A/j, (but 

small compared to Ep and /csT ). We have calculated 
the demagnetization with finite spin relaxation time of 
the holes using the properly down-scaled transition rates 
from Eq. ([27]) with finite A/i. The results are shown in 
Figure [5] The effect of spin bottleneck on demagnetiza- 
tion is clearly very weak for r sr <10 fs. For t si =100 fs the 
total demagnetization is diminished by ^50%. These re- 
sults agree with estimates from the inequality (|38| , which 
is fulfilled for 5^1 meV and r sr =10 fs. 

In the presence of the strong spin-orbit interaction, the 
hole spin relaxation is expected to be very fast, occurring 
on the momentum scattering time-scaled For example, 
in pure GaAs the spin relaxation of holes was measured to 
be about 100 fsj£L The calculated^ momentum scatter- 
ing of holes in GaMnAs at low temperature is about 10 fs, 
and it is expected to be shorter when the holes are highly 
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Markovian regime, 78 and r sr is a typical time-scale of the 
spin-dependent correlated relaxation dynamics. 

In Fig. [8] we can see that for r sr <100 the spin bottle- 
neck effect is qualitatively unimportant, and it becomes 
negligible for r sr =10 fs. Without the bottleneck, we ob- 
tain a typical magnitude of demagnetization in (III,Mn)V 
within a picosecond to be of the order of 10%. For the 
same carrier density, initial T c =1000 K and te=0.5 ps, 
the drop in magnetization is more pronounced in GaM- 
nAs (20% demagnetization) than in InMnAs (10%). A 
discussion of the connection between these calculations 
and experiments follows below. 



COMPARISON WITH EXPERIMENTAL 
RESULTS 



FIG. 8: (Color online) Demagnetization of the Mn spin in (a) 
GaMnAs and (b) InMnAs. The solid lines are calculated using 
the exact transition rates at Afi=0 (dynamical spin polariza- 
tion of carriers is neglected). The dashed lines come from the 
approximate calculation including the finite hole spin relax- 
ation time r sr (see text for details). The initial temperature 
of holes if T c = 1000, the hole density is p=3 x 10 20 cm -3 , the 
exchange integral /3=— 50 meV-nm 3 , and the molar fraction 
of Mn is 1=0.054 (initial spin splitting A=-0.15 eV). The 
energy relaxation time te is 0.5 ps. 

excited. There is one caveat in the case of exchange spin- 
split bands in (III,Mn)V. The lifting of the heavy-light 
hole degeneracy by confinement in quantum wells makes 
the spin relaxation time longer , 73 ' 74 with measured val- 
ues of 4 ps for holes close to T point in modulation-doped 
GaAs quantum wells. 75 An analogous effect is expected 
in exchange-split bands for holes with small wave- vectors. 
However, for the carrier densities and temperatures under 
consideration we are mostly interested in spin relaxation 
of holes with quite large wavevectors (|k|«l— 3 urn -1 ). 
For such large |k| the effect of lifted degeneracy on spin 
relaxation becomes much weaker^ For example, at a 
wavevector of this magnitude in the direction perpendic- 
ular to the magnetization the band-mixing terms in the 
Luttinger Hamiltonian overcome the p-d interaction, and 
the spin splitting becomes negligible. Taking all this into 
account it is reasonable assume that in an excited and 
disordered sample the hole spin relaxation time should 
be shorter than in pure GaAs. For r sr <100 fs we can see 
in Fig. [8] that the spin bottleneck effect is qualitatively 
unimportant, and it becomes negligible for r sr = 10 fs. 

In the experiments on (III,Mn)V spin lifetimes of elec- 
trons between 1 ps in highly excited InMnAs 7 -^ and 30 ps 
in much more weakly excited GaMnAs 7 ^ were seen. To 
the best of our knowledge, no signal attributable to hole 
spin relaxation has ever been seen on a time-scale of the 
temporal resolution of the experiments. If r sr <10 fs, the 
measurements would be particularly challenging due to 
time-resolution constraints on the laser pulse time. Also, 
at such time-scale the carrier dynamics is in the non- 



We can estimate the initial temperature of carriers 
in experiments on InMnAs from Ref. [Ill in the follow- 
ing way. There are two time-scales in the demagnetiza- 
tion: an ultrafast one (<1 ps) and a subsequent much 
slower (~100 ps) demagnetization, due to spin-lattice re- 
laxation. For example, for pump flucncc of 3 mJ/cm 2 , 
an ultrafast quenching of ~ 50% of magnetization was 
followed by a complete demagnetization 100 ps later. On 
the latter time-scale we can safely use a thermal descrip- 
tion of the spin system. The full demagnetization means 
that the localized spin temperature T s had risen above 
the Curie temperature, which was 50 K in the sample 
used. This tells us that 100 ps after the excitation the 
lattice temperature T; was at least 50 K. The energy (per 
atom) , which had to be transferred into the lattice to heat 
it up from initial temperature Ti (before the excitation) 
up to the final T; is given by 

AE l= / c L {T)dT , (39) 

JTi 

where c; is the specific heat of the lattice (per atom), 
which in the temperature range of interest is given by 
the Debye formula: 

ct(T) = — k B (-) , (40) 

where 0^=280 K is the Debye temperature of InAs. 

The energy transferred into the lattice needs now to 
be related to the energy deposited initially by the light 
pulse into the valence band. In the process of absorption, 
the energy of the pump photon (hu=0.6 eV) is split into 
a kinetic energy of a created hole and an energy of an ex- 
cited state within a gap (we disregard the photoelectrons 
because the low density of states in the conduction band 
makes this excitation channel insufficient to explain the 
experiments) . The Fermi level before the excitation is sit- 
uated between 0.1 and 0.3 eV below the top of the valence 
band for hole densities of 10 20 and 3 x 10 20 cm -3 , respec- 
tively. Thus, the energies of the newly created holes are a 
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sizable fraction of the photon energy. As a rather safe es- 
timate we will take the fraction of the pump pulse energy 
imparted to the holes as i?=l/4. 

Now we can derive an estimate of the carrier temper- 
ature T c (0) just after the absorption of the pump pulse. 
We assume that after 100 ps all the absorbed energy has 
been transferred to the lattice. The observed spin tem- 
perature T s is larger than 50 K and gives a lower bound 
for the final lattice temperature T;. We calculate the en- 
ergy given to the lattice using Eqs. ([39]) and (|40|) . with 
initial T;=10 K and final T;=50 K. A fraction R of this 
energy was an excess kinetic energy of holes in the valence 
band after the excitation. Using the density of states of 
InMnAs with A=0.15 eV calculated before, we obtain the 
temperature of carriers which gives such excess energy. 
In this way we obtain an estimate of the initial T e , which 
we find to be between 1500 and 1000 K for hole densities 
changing between 10 20 and 3 x 10 20 cm 3 . The analogous 
calculation in GaMnAs gives T c smaller by a factor of 2 
(due to the larger density of states). These results jus- 
tify our use of typical initial To =1000 K throughout this 
article. 

Let us now address the corresponding demagnetization 
in InMnAs for pump flucnce of 3 mJ/cm 2 . The hole con- 
centration before the excitation is about 3 x 10 20 cm -3 , 
estimated from critical temperature of 50 K and magne- 
tization measurements After the excitation the total 
number of holes is larger but, because the initial p is not 
certain, we will simply use a value of total p=3 x 10 20 
cm" 3 . The measured sample was a 25 nm thick layer of 
InMnAs grown on GaSb. For such thickness, the confine- 
ment leads to formation of hole subbands, with typical 
energy spacing between them of the order of 10 meV. 
This confinement energy is much less than the disorder 
broadening and the thermal spread of photoexcited car- 
riers (both ~0.1 eV), so that our calculation of spin- flip 
transition rates using the bulk band-structure should be 
a good approximation. In the previous section, we have 
seen that for such p=3 x 10 20 cm -3 the sub-picosecond 
drop in magnetization is about 10%. Experimentally a 
50% drop in Kerr signal was observed within ^200 fs after 
the pump pulse. The step-like character of the magneti- 
zation drop is not reproduced by our theory, which pre- 
dicts a smoother demagnetization. However, the carrier- 
induced artifacts (termed "dichroic bleaching" in Ref. 0) 
obscuring the magnetization dependence of the magne- 
tooptical signal are possible at very short time-scale. The 
magnitude of the total demagnetization is in qualitative 
agreement with our calculations. 

For pump fluence of 10 mJ/cm 2 , for which a complete 
quenching of ferromagnetic order is observed, the role of 
multi-photon absorption processes becomes pronounced, 
and photocarriers arc created in a very large region of the 
Brillouin Zone (as in experiments from Ref. 79, where 
a pump with Hlo—3.1 eV was used in GaMnAs). Our 
kp model with "rigid" spin splitting (i.e. k-independent 
exchange constant) is not applicable far away from the T 
point. Also, for very high excitation the assumption of 



quasi-thermal distribution of holes can fail. 

In GaMnAs the experiments show a similar behavior- 
after 1 picosecond, with 30% magnetization drop for flu- 
ence of ~8 mJ/cm 2 . This value is very close to the pre- 
dictions of our theory for p~3 x 10 20 cm~ 3 . For the same 
fluence InMnAs was already nearly completely demag- 
netized. Although our theory predicts larger demagne- 
tization in GaMnAs than in InMnAs, in order to make 
meaningful comparisons one has to achieve comparable 
excitation parameters (total p and T D after the pulse). 
Also, the efficiency of excitation of holes can be different 
in the two materials. 



VI. SUMMARY 

We have presented an investigation of the ultrafast de- 
magnetization induced by an incoherent light excitation 
in a system described by the sp-d model. The physical 
picture of demagnetization is very transparent, with sp-d 
interaction providing a mechanism of spin transfer from 
the localized spins into the system of excited carriers. 
This process is closely related to effects known from sys- 
tems of electronic and nuclear spins coupled by hyperfine 
contact interaction, and thus we have termed it the in- 
verse Overhauser effect. The demagnetization basically 
comes from very fast T\ relaxation of the localized spins 
due to sp-d interaction with the hot carriers' spins. In a 
simple one-band model the rate of spin flip scattering of 
itinerant and localized spins is proportional to the square 
of the sp-d exchange constant, the product of electronic 
densities of states at Fermi energy for two spin directions, 
and the temperature of the carriers. Due to the carrier 
temperature dependence the rapid demagnetization can 
occur only for a time of the order of energy relaxation 
time te- Thus, the sub-picosecond time scale of the de- 
magnetization process simply comes from the character- 
istic time of carrier-phonon interaction. The total magni- 
tude of demagnetization during the te time depends only 
on the spin-flip scattering efficiency if the carrier system 
is a good spin sink, i.e. the carrier spin relaxation time 
r sr is shorter than the time-scale on which the dynamic 
spin polarization of carriers builds up. A long T sr leads 
to spin bottleneck that suppresses the demagnetization. 
This general theory can be used as an admittedly rough 
model for demagnetization process in transition metals 
(treated within s-d model), as it captures all the crucial 
aspects of the problem^ 

A large part of the discussion was aimed at the specific 
cases of (III,Mn)V semiconductors. There, a non-trivial 
band structure with spin-orbit interaction has to be taken 
into account. We have performed the calculations of spin- 
flip transition rate using an effective k-p Hamiltonian ap- 
proach. An approximate calculation of the effects of dy- 
namical spin polarization of holes has been used to argue 
that the holes can be treated as a perfect sink for the an- 
gular momentum transferred from the localized spins. A 
qualitative agreement with demagnetization experiments 
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in these materials was obtained: the theory shows that 
the magnetization can drop by ~10% within the energy 
relaxation time of the holes. More theoretical and ex- 
perimental work on electronic and optical properties of 
(III,Mn)Vs will be necessary in order to better reconcile 
experiments and theory. 
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